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Abstract. Let R denote dyadic rectangles in the unit cube [0, 1]^ in three dimensions. Let Hr 
be the L°°-normali2ed Haar function whose support is R. We show that for all integers n > 1 
and choices of coefPicients qr € {±1}, we have 



Y aRh.R 



> TL 



9/8 



|R|=2 — 
|R, |>2-"''^ 

The trivial lower bound is n, and the sharp lower bound would be n^^^. This is the best expo- 
nent known to the authors. This inequality is motivated by new results on the star-Discrepancy 
function in all dimensions d > 3. 



1. Introduction 

We are motivated by the classical question of irregularities of distribution [2] and recent results 
which give new lower bounds on the star-Discrepancy in all dimensions d > 3 ]. We recall 
these results. 

Given integer N, and selection V of N points in the unit cube [0, 1]'^, we define a Discrepancy 
Function associated to V as follows. At any point x G [0, 1]'^, set 

Dn(x) =jj(Pn[0,x))-N|[0,x)|. 

Here, by [0, x) we mean the d-dimensional rectangle with left-hand corner at the origin, and 
right-hand corner at x G [0, 1]'^. Thus, if we write x = (xi , . . . , Xd) we then have 

d 

j=i 

At point X we are taking the difference between the actual and the expected number of points in 
the rectangle. Traditionally, the dependence of Dm on the selection of points V is only indicated 
through the number of points in the collection V. We mention only the main points of the subject 
here, and leave the (interesting) history of the subject to references such as [ ]. 

The result of Klaus Roth [ ] gives a definitive average case lower bound on the Discrepancy 
function. 
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K. Roth's Theorem. For any dimension d > 2, we liave tlie following estimate 

||Dn||2> (logN)('i-^'/^ 

The same lower bound holds in all L^^, 1 < p < oo, as observed by Schmidt [8]. But, the L°° 
infinity estimate is much harder. In dimension d — 2 the definitive result was obtained by Schmidt 
again [9]. 

Schmidt's Theorem. We have the estimates below, valid for all collections V C [0, 1]^; 

IIDnIIoo > logN. 

The L°° estimates are referred to as star-Discrepancy bounds. Extending and greatly simplifying 
an intricate estimate of Jozef Beck ['], some of these authors have obtained a partial extension 
of Schmidt's result to all dimensions d > 3. 

Theorem 1.1. /[4, 5]/ For dimensions d > 3 there is an r| = ri(d) > for which we have the 
inequality 

||DN|U>(logN)('i-^'/^+^ 
That is, there is an r\ improvement in the Roth exponent. 

As explained in these references, the analysis of the star-Discrepancy function is closely related 
to other questions in probability theory, approximation theory, and harmonic analysis. We turn to 
one of these, the simplest to state question, which is central to all of these issues. We begin with 
the definition of the Haar functions. 

In one dimension, the dyadic intervals of the real line R are given by 

P = {[j2\ + 1)2^) : j,kGZ}. 

Any interval I is a union of its left and right halves, denoted by I|eft/right. which are also dyadic. 
The Haar function Hi associated to I, or simply Haar function is 

Hi = -llleft + llright 

Note that for dyadic intervals J £ I, the Haar function hj is completely supported on a set where 
hi is constant. This basic property leads to far-reaching implications that we will exploit in these 
notes. 

In higher dimensions d > 2, we take the dyadic rectangles to be the tensor product of dyadic 
intervals in dimension d: 

^ ^ X . .. X : Ri,...,Rd G P} . 
The Haar function associated to R G Pa is likewise defined as 

d 

Hr(xi, . . .,Xd] = J^HrjIxj] , R = Ri X • ■ ■ X Rd. 
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While making these definitions on all of R'^, we are mainly interested in local questions, thus 
rectangles R C [0, 1]'^ are always dyadic rectangles R G V^. Namely, we are mainly interested in 
the following conjectural reverse triangle inequality for sums of Haar functions on L°°: 

Conjecture: The Small Ball Inequality. For dimensions d > 3, there is a constant Ca so that 
for all integers n > 1 , and constants {a^ : |R| = 2^^ , R C [0, 1]'^}, we have 



(1.2) 



n 



(d-2)/2 



R|>2-'^ 
RC[0,1]'* 



> Ch2- 



|R|=2-" 
Rc[0,l]'i 



I^R 



We are stating this inequality in its strongest possible form. On the left, the sum goes over all 
rectangles with volume at least 2^^, while on the right, we only sum over rectangles with volume 
equal to 2^^. Given the primitive state of our knowledge of this conjecture, we will not insist on 
this distinction below. 

For the case of d = 2, (1.2) holds, and is a Theorem of Talagrand [10]. (Also see [6, 8, 11]). 

The special case of the Small Ball Inequality when all the coefficients Qr are equal to either 
— 1 or +1 we refer to as the 'Signed Small Ball Inequality.' Before stating this conjecture, let us 
note that we have the following (trivial) variant of Roth's Theorem in the Signed case: 



Y <lR ■ Hr 



|R|=2-'^ 
Rc[0,l]d 



> Tl 



(d-l)/2 



aR G{±1}, 



The reader can verify this by noting that the left-hand side can be written as about n'^^^ orthogonal 
functions, by partitioning the unit cube into homothetic copies of dyadic rectangles of a fixed 
volume. The Signed Small Ball Inequality asserts a 'square root of n' gain over this average case 
estimate. 

Conjecture: The Signed Small Ball Inequality. For coefficients Qr G {±1}, 



Y ■ It-R 



|Rh2— 
RC[0,1]'' 



Here, is a constant that only depends upon dimension. 

We should emphasize that random selection of the coefficients shows that the power on n on 
the right is sharp. Unfortunately, random coefficients are very far from the 'hard instances' of the 
inequality, so do not indicate a proof of the conjecture. 

The Signed Small Ball Conjecture should be easier, but even this special case eludes us. To 
illustrate the difficulty in this question, note that in dimension d = 2, each point x in the unit 
square is in rL+ 1 distinct dyadic rectangles of volume 2^^. Thus, it suffices to find a single point 
where all the Haar functions have the same sign. This we will do explicitly in § 2 below. 

Passing to three dimensions reveals a much harder problem. Each point x in the unit cube is 
in about rectangles of volume 2^^, but in general we can only achieve a jv'^^ supremum norm. 
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Thus, the task is to find a single point x where the number of pluses is more than the number of 
minuses by n?^^. In percentage terms this represents only a nr^^^-percent imbalance over equal 
distribution of signs. 

The main Theorem of this note is Theorem 4.1 below, which gives the best exponent we are 
aware of in the Signed Small Ball Inequality. The method of proof is also the simplest we are 
aware of. (In particular, it gives a better result than the more complicated argument in [3]). 
Perhaps this argument can inspire further progress on this intriguing and challenging question. 

The authors thank the anonymous referee whose attention to detail has helped greater clarity 
in our arguments. 

Dedication to Walter Philipp. One of us was a PhD student of Walter Philipp, the last of seven 
students. Walter was very fond of the subject of this note, though the insights he would have into 
the recent developments are lost to us. As a scientist, he held himself to high standards in all his 
areas of study. As a friend, he was faithful, loyal, and took great pleasure in renewing contacts 
and friendship. He is very much missed. 

2. The Two Dimensional Case 

This next definition is due to Schmidt, refining a definition of Roth. Let f G N'^ be a partition 
of n, thus f = (ti , . . . , Tii), where the tj are non negative integers and |f| := ^^=1 = "n.- 
Denote all such vectors at M.^. ('H' for 'hyperbolic.') For vector f, let T^f be all dyadic rectangles 
R such that for each coordinate 1 < t < d, |Rt| — 2^'^^ . 

Definition 2.1. We call a function f an r-function with parameter r if 

Y. ^K^^' e{±l}. 

We will use to denote a generic r-function. A fact used without further comment is that fi=]. 

Note that in the Signed Small Ball Inequality, one is seeking lower bounds on sums Y.\T\=n'^?- 
There is a trivial proof of the two dimensional Small Ball Inequality. 

Proposition 2.2. Ttie random variables f(j,n-j), < ] < n are independent. 

Proof. The sigma-field generated by the functions {f(k,n-k) '■ < k < j} consists of dyadic 
rectangles S = Si x Si with |Si| = 2^' and IS2I = 2^^. On each line segment Si x {X2}, f(j,n-)) 
takes the values ±1 in equal measure, so the proof is finished. 

□ 

We then have 
Proposition 2.3. In the case of two dimensions, 

k=0 
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Proof. Note that 

' f (k,n-k] = n + 1 ] = P(f (k,n-k) = 1 VO < k < n) = 2-"-^ 

' k=0 

□ 

It is our goal to give a caricature of this arugment in three dimensions. See §5 for a discussion. 

3. Elementary Lemmas 

We recall some elementary Lemmas that we will need in our three dimensional proof. 

Paley-Zygmund Inequality. Suppose that Z is a positive random variable witti EZ = \L), 
EZ^ = \il Then, 

(3.1) P(Z>^i/2)>l4. 
Proof. 

= EZ = EZlz<^,/2 + EZlz>^,/2 
< m/2+^2P(Z> 

Now solve for P(Z > |J.i/2). □ 

Second Paley-Zygmund Inequality. For all pi > 1 there is a Qi > ^ so that for all random 
variables Z which satisfy 

(3.2) EZ = 0, ||Z||2<l|Zl|4<Pil|Z||2 
we have the inequality ¥[Z > P2IIZII2) > P2- 

Proof Let Z+ := Zlz>o and Z := -Zlz<o, so that Z = Z+ - Z . Note that EZ = forces 
EZ+ =EZ_. And, 

al ■- EZ^ = EZ^ + EZ^ , 
:= EZ'* = EZl + EZi . 

Suppose that the conclusion is not true. Namely P(Z > P2cr2) < P2 for a very small P2. it 
follows that 

EZ+ < EZ+lz+<p2a2 +EZ+lz+>p2a2 

< P20-2 + P(Z > P20-2)^/'0-2 < 2py'cT2 , 

1 /2 

for P2 < 1 . He nee EZ_ — EZ+ < 2p2 (Ji- It is this condition that we will contradict below. 
We also have 

EZ^ < EZ^+l{z,<p,a,} +EZil{z,>p,a,} 

/ 2 2 , 1/2 7 
< P20-2 + P2 O4 
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So for P2 < (4pi)^'*, we have EZ^ < ^ct|. 

It follows that we have EZi > jGj, and EZi < Picr^. So by (3.1), we have 

P(Z_ > paai) > P3 

where ps is only a function of pi. But this contradicts EZ_ < 2p2'^^ff2, for small p2, so finishes 
our proof. □ 

The Paley-Zygmund inequalities require a higher moment, and in application we find it conve- 
nient to use the Littlewood-Paley inequalities to control this higher moment. Let J-o, J-"] , . . . , J-j a 
sequence of increasing sigma-fields generated by dyadic intervals, and let dt, 1 < t < T be a mar- 
tingale difference sequence, namely E(dt : J^t-i) = for all t = 1 , 2, . . . , T. Set f = dt. 
The martingale square function of f is S(f)^ := ^{=1 '^t- The instance of the Littlewood-Paley 
inequalities we need are: 

Lemma 3.3. With the notation above, suppose that we have in addition that the distribution of 
dt is conditionally symmetric given J^t-i ■ By this we mean that on each atom A of J-'t-i , the the 
distribution of dtlA is equal to that of — dtlA- Then, we have 

(3.4) ||f||4^||S(f)||4. 

Proof. The case of the Littlewood-Paley for even integers can be proved by expansion of the 
integral, an argument that goes back many decades, and our assumption of being conditionally 
symmetric is added just to simplify this proof. Thus, 

l|f||^= Y_ ^tldu- 

l<tl,t2,t3,t4<T U=l 

We claim that unless the integers 1 < ti,t2,t3,t4 < T occur in pairs of equal integers, the 
expectation on the right above is zero. This claim shows that 

l<tl,t2<T 

It is easy to see that this proves the Lemma, namely we would have 

iis(f)iii<iifii^<fiis(f)iir 

Let us suppose ti < t2 < t3 < 14. If we have strictly less than t4, then 

4 3 
Efl dt. = Ef] dt, ■ E(dt, : Tt,) = . 

u=l u=l 

If we have ti < t2 = t3 = 14, then by conditional symmetry, E(d^^ : J-'t, ) = 0, and so we have 

4 

Ej^dt. =Edt, - Eidl : -Ft,) =0. 

u=l 
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If we have ti < t2 < ts = t4, the conditional symmetry again implies that E(d.t2 ' '■ J^u) — 0, 
so that 

4 

e]^ dt, = Edt, ■ E(dt, ■ d^3 : j;, ) = . 

Thus, the claim is proved. 



□ 



We finish this section with an elementary, slightly technical. Lemma. 



Lemma 3.5. Let Fq^ J-"] , . . . , J-'q a sequence of increasing sigma-fields. Let Ai , . . . , Aq be events, 
with At G J-t. Assume ttiat for some < y < 1 , 

(3.6) e(1a, : J't-^) >y, 1 < t< q 
We then have that 

(3.7) p(nAt)>y'>. 
More generally, assume that 

(3.8) p(U{e(1a. : ^t-i) <y})<ly'. 



t=l 



Then, 



(3.9) p(nAt) > ^y^ 

Proof To prove (3.7), note that by assumption (3.6), and backwards induction we have 

q X q 
HAt ^eHIa, 

t=i ^ t=i 
q-i 

= e]^1a, xE(1a^ : J-q_i 
t=i 
q-i 

>yE]^lA, 

t=i 



> y^ . 

To prove (3.9), let us consider an alternate sequence of events. Define 

|3t := {E(At : J-t-i) <y}. 
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These are the 'bad' events. Now define At := At U (3f By construction, the sets At satisfy (3.6). 
Hence, we have by (3.7), 



But, now note that by (3.8), 



nAt)=p(nAt)-p(upt 

t=l ^ H=l ^ H=l 

>y''-;-y'^>;-y^ 



11-1 

□ 



4. Conditional Expectation Approach in Three Dimensions 

This is the main result of this note. 
Theorem 4.1. For laRl = 1 for a\\ R, we have the estimate 



|Rh2— 
|R, |>2-'^/2 



We restrict the sum to those dyadic rectangles whose first side has the lower bound |Ri | > 2 

Heuristics for our proof are given in the next section. The restriction on the first side lengths of 
the rectangles is natural from the point of view of our proof, in which the first coordinate plays a 
distinguished role. Namely, if we hold the first side length fixed, we want the corresponding sum 
over R to be suitably generic. Let 1 -C q <^ n be integers. The integer q will be taken to be 
q ~ TL^/'*. Our 'gain over average case' estimate will be y/q ~ m}^^. While this is a long way 
from n'^^, it is much better than the explicit gain of 1 /24 in [3]. 

We begin the proof. Let J-'t be the sigma field generated by dyadic intervals in [0,1] with 
|I| = for 1 < t < ^q. Let It := {f : (t - l)n/q < r^ < tn/q}. Note that the size 

#It ~ Tt^/q. Let ff be the r-functions specified by the choice of signs in Theorem 4.1. Here is a 
basic observation. 

Proposition 4.2. Let I G J^t- The distribution of {ff : f G It} restricted to the set I x [0, 1]^ 
with normalized Lebesgue measure is that of 

{fs- : |s| =n - [tn/qj , < s, < n/q}, 

where the f ^ are some r-functions. The exact specification of this collection depends upon the 
atom in T^. 

Proof. An atom I of Tt are dyadic intervals of length l^LWqJ pgr f ^ restricted to 

I X [0, 1]^, with normalized measure, is an r-function with index 

(ri - [tn/q J, r2,r3) . 

The statement holds jointly in rG It so finishes the proof. □ 
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Define sum of 'blocks' of ff as 



(4.3) 



Tgit 

nt:= ff-fs 



T,=S, 



The sums fit P'^Y ^ distinguished role in our analysis, as revealed by the basic computation of 



a square function in (4.8) and the fundamental Lemma 4.9. Let us set 
< t < q/2. 

We want to show that for q as big as cn^/^, we have 



|Bt||2 — TL^/q, for 



(4.4) 

In fact, we will show 



q/2 
t=i 



q/2 



P n{Bt>n/v/q} >0, 



t=i 



from which (4.4) follows immediately. 

Note that the event |Bt > n/ y^} simply requires that Bt be of typical size, and positive, 
that is this event will have a large probability. Clearly, we should try to show that these events are 
in some sense independent, in which case the lower bound in (4.4) will be of the form e^'-'', for 
some C > 0. Exact independence, as we had in the two-dimensional case, is too much to hope 
for. Instead, we will aim for some conditional independence, as expressed in Lemma 3.5. 

There is a crucial relationship between Bt and \~]^, which is expressed through the martingale 
square function of Bt, computed in the first coordinate. Namely, define 

2 

(4.5) S(Bt)^ ■= II ^ 

jejt f : T,=j 

: (t — l)n/q < s < tn/q}. 



where Jt = {s G I 
Proposition 4.6 

(4.7) 
(4.8) 



We have 

s(Bt]^ = o^ + nt, 

S(Bt : ^t) = oi + E(nt : -^t) • 
By construction, we have ''it ^ Ti^/q, for < t < q. 
Proof. In (4.5), one expands the square on the right hand side. Notice that this shows that 



S(Bt 



|rh|s|=n 
ri=si gjt 



10 D. BILYK, M.T. LACEY, I. PARISSIS, AND A. VAGHARSHAKYAN 

We can have f = s for ''It choices of f. Otherwise, we have a term that contributes to fit- The 
conditional expectation conclusion follows from (4.7) □ 

The next fact is the critical observation in [3-5] concerning coincidences, assures us that 
typically on the right in (4.7), that the first term ~ n^/q is much larger than the second fit- 
See [ , 4.1, and the discussion afterwords]. 

Lemma 4.9. We have the uniform estimate 

Here, we are using standard notation for an exponential Orlicz space. 

Remark 4.10. A variant of Lemma 4.9 holds in higher dimensions, which permits an extension of 
Theorem 4.1 to higher dimensions. We return to this in §5. 

Let us quantify the relationship between these two observations and our task of proving (4.4). 

Proposition 4.11. There is a universal constant t > so that defining the event 

(4.12) rt:={E(n^ : T,y/'<xn'/q] 
we have the estimate 

(4.13) P(Bt>T-n/v^ : > Tip, . 

The point of this estimate is that the events Ft will be overwhelmingly likely for q <C n. 

Proof This is a consequence of the Paley-Zygmund Inequalities, Proposition 4.2, Littlewood- 
Paley inequalities, and (4.8). Namely, by Proposition 4.2, we have E(Bt : J-'t) = 0, and the 
conditional distribution of Bt given J^t is symmetric. By (4.8), we have 

E(B^ : J-t)=S(Bt : T,) = al + E{nt : -^t) • 
We apply the Littlewood-Paley inequalities (3.4) to see that 

E(B^ : J-J <E(S(Bt : J^t)^ : J^t) 

= a4 + 2a?E(nt : -^t)+E(nt : -^t) • 

The event Ft gives an upper bound on the terms involving fit above. This permits us to 
estimate, as Ft G J^t. 

E(B^ : Ft)^/^-atlr,| <Tn/v^, 
but at — tl/ ^/ai, so we have E(Bt : Ft)^/^ ~ n/ y/q. Similarly, 

E(B^ : Ft)^/^<(Tt + (Ty'|E(nt : J't]\'^' + E[ni : J'tY^' 

< (1 +T)cTt. 

Hence, we can apply the Paley-Zygmund inequality (3.2) to conclude the Proposition. □ 
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By way of explaining the next steps, let us observe the following. If we have 
(4.14) E(lr, : j;) > t: a.s. (xi), 1<t<q/2, 

then (4.13) holds, namely P(Bt(-, xi, X3) > t ■ n/ : Ft) > t almost surely. Applying 
Lemma 3.5, and in particular (3.7), we then have 



"1 



q/2 

n{Bt(-,X2,X3]>Tn/v/q})>T^/^ 
t=i 



Of course there is no reason that such a pair (xijXa) exits. Still, the second half of Lemma 3.5 
will apply if we can demonstrate that we can choose X2,X3 so that (4.14) holds except on a set, 
in the Xi variable, of sufficiently small probability. 

Keeping (3.8) in mind, let us identify an exceptional set. Use the sets Ft as given in (4.12) to 
define 



(4.15) 



E := 





"q/2 


|(X2,X3) : p., 


_t=i 



> 



exp(^-Ci(n/q) 



1/3 



Here, Ci > will be a sufficiently small constant, independent of n. Let us give an upper bound 
on this set. 



■ "2 



,X3(E) < exp(ci(rL/q) 



\l/3 



q/2 
t=i 



q/2 



t=i 

qexp(ci(Ti/q)^/3) ■exp|^-[T(nVq)||E(nt : -^t 
< q exp([ci - C2T^/^) ■ (n/q)^/^ 



< 



1/2, 



-1 

lcxp(L2/3) 



2/3' 



(4.16) 

Here, we have used Chebyscheff inequality. And, more importantly, the convexity of conditional 
expectation and L^-norms to estimate 



1/2, 



lcxp(L2/3) 



by Lemma 4.9. The implied constant is absolute, and determines the constant C2 in (4.16). For 
an absolute choice of Ci, and constant t', we see that we have 



■ "2 



,,3(E) <exp(-T'(n/q) 



l/3l 



We only need Px2,x3(t) < j, but an exponential estimate of this type is to be expected. 
Our last essential estimate is 
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Lemma 4.17. For < k < 1 sufficiently small, q < Kn^^"^, and (x2,X3) ^ E, we have 



(n{Bt(-,X2,X3) >Tn/v/q}) > 



Assuming this Lemma, we can select (x2,X3) ^ E. Thus, we see that there is some (xi,X2,X3) 
so that for all 1 < t < q/2 we have Bt(xi,X2,X3) > ttl/ ^/q, whence 

q/2 

^Bt(xi,X2,X3] >j-n^. 
t=i 

That is, (4.4) holds. And we can make the last expression as big as > n^/^. 

Proof If (x2,X3) ^ E, bring together the definition of E in (4.15), Proposition 4.11, and Lemma 3.5. 
We see that (3.9) holds (with y = t, and the q in (3.9) equal to the current q/2) provided 

>exp(-ci(n/q)^/^) . 

But this is true by inspection, for q < Kn^/"*. 

□ 



5. Heuristics 

In two dimensions. Proposition 2.3 clearly reveals an underlying exponential-square distribution 
governing the Small Ball Inequality. The average case estimate is n^^^, and the set on which the 
sum is about n (a square root gain over the average case) is exponential in n. 

Let us take it for granted that the same phenomena should hold in three dimensions. Namely, 
in three dimensions the average case estimate for a signed small ball sum is n, then the event 
that the sum exceeds n^''^ (a square root gain over the average case) is also exponential in n. 
How could this be proved? Let us write 

n/2 

|R|=2-'^ |f|=n j=0 

|R, |>2-" ri<n/2 

\r\=n 

Here we have imposed the same restriction on the first coordinate as we did in Theorem 4.1. 
With this restriction, note that each |3j is a two-dimensional sum, hence by Proposition 2.2, a 
sum of bounded independent random variables. It follows that we have by the usual Central Limit 
Theorem, 
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for a fixed constant c. If one could argue for some sort of independence of the events {|3j > c y/n} 
one could then write 

P(H > cn^/^) > P f|{(3j > c v^} > , 

for some e > 0. This matches the 'exponential in n' heuristic. We cannot implement this proof 
for the (3j, but can in the more restrictive 'block sums' used above. 

We comment on extensions of Theorem 4.1 to higher dimensions. Namely, the methods of this 
paper will prove 

Theorem 5.1. For [qrI = 1 for all R, we have the estimate estimate in dimensions d > 4; 

,(d-l)/2+l/4d 



|R|=2-'^ 
|Rl|>2-'^/2 



We restrict the sum to those dyadic rectangles whose first side has the lower bound |R| > 2 

This estimate, when specialized to d = 3 is worse than that of Theorem 4.1 due to the fact that 
the full extension of the critical estimate Lemma 4.9 is not known to hold in dimensions d > 4. 
Instead, this estimate is known. Fix the coefficients aR G {±1} as in Theorem 5.1, and let ff be 
the corresponding r-functions. For 1 -C q -C n, define It as above, namely {f : \r\ — n , vi E Jt}. 
Define Ht in (4.3). The analog of Lemma 4.9 in dimensions d > 4 are 

Lemma 5.2. In dimensions d> 4 we have the estimate 

l|ntllexp(L^/.--i.)<n(^'^-31/Vv^. 

See [4, Section 5, especially (5.3)], which proves the estimate above for the case of q = 1. 
The details of the proof of Theorem 5.1 are omitted, since the Theorem is at this moment only 
a curiosity. 

It would be quite interesting to extend Theorem 5.1 to the case where, say, one-half of the 
coefficients are permitted to be zero. This result would have implications for Kolmogorov entropy 
of certain Sobolev spaces; as well this case is much more indicative of the case of general coeffi- 
cients qr. As far as the authors are aware, there is no straight forward extension of this argument 
to the case of even a small percentage of the Qr being zero. 
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